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ABSTRACT 



Context. The potential role of magnetic fields and cosmic ray propagation for feedback processes in the early Universe can be probed 
by studies of local starburst counterparts with an equivalent star-formation rate. 

Aims. In order to study the cosmic ray propagation, determine the magnetic field strength and the dominant loss processes in the 
nearby prototypical starbursting galaxy M82 a multi-frequency analysis at four radio wavelengths is presented. 
Methods. Archival data from the WSRT was reduced and a new calibration technique introduced to reach the high dynamic ranges 
needed for the complex source morphology. This data was combined with archival VLA data, yielding total power maps at A3 cm, 
A6 cm, A22 cm and A92 cm. 

Results. The data shows a confinement of the emission at wavelengths of A3/A6 cm to the core region and a largely extended halo 
reaching up to 4kpc away from the galaxy midplane at wavelengths of A22/A92cm up to a sensitivity limit of 90 \ijy and 1.8 mJy 
respectively indicating different physical processes in the core and halo regions. The results are used to calculate the magnetic field 
strength in the core region to 98 \iG and to 24 u.G in the halo regions. From the observation of ionisation losses the filling factor of 
the ionised medium could be estimated to 2%. This leads to a revised view of the magnetic field distribution in the core region and 
the propagation processes from the core into the halo regions. 

Conclusions. We find that the radio emission from the core region is dominated by very dense HH-regions and supernova remnants, 
while the surrounding medium is filled with hot X-ray and neutral gas. Cosmic rays radiating at frequencies higher than 1.4 GHz are 
suffering from high synchrotron and inverse Compton losses in the core region and are not able to reach the halo. Even the cosmic 
rays radiating at longer wavelengths are only able to build up the observed kpc sized halo, when several starbursting periods are 
assumed where the photon field density varies by an order of magnitude. These findings together with the strong correlation between 
Ha, PAH + , and our radio continuum data suggests a magnetic field which is frozen into the ionised medium and driven out of the 
galaxy kinematically. 

Key words. Galaxies: individual: M82 - Galaxies: starburst - Galaxies: magnetic fields - Galaxies: halos - Techniques: interferometric 



1. Introduction 

M82 is a nearby prototypical starburst galaxy with a distance 
of ~3.5Mpc ^Jacobs et al.|20 09 ) in the M81 group of galaxies. 
It hosts a complicated distribution of different gas phases in its 
halo as well as in its violently starbursting core. Material from 
the core region is accelerated into the halo via a galactic super- 
wind as well as dragged out of the galaxy by tidal interaction. A 
prominent large scale biconical outflow has been observed in Ha 
I tiDevine & Bally|1999[|Ohyama et al.|2002| ) as well as in X-rays 
( Strickland & Heckman 2009 ) Dust, neutral gas phases and even 
molecular material are mixed in the ionised phase and seem to 
be swept along by this superwind (e.g. CO: [Walt er et al.| (|2002|i j 



H 2 : Veilleux et aL] (2009 ) ; Dust:|Ohyama et al. (2002i; Leeuw 



& Robson (2009)). Large-scale molecular and atomic streamers 
are anchored at the edges of the galactic disk and therefore show 
a completely different distribution (Wa lter et al.|2002 Yun et al. 
1993b) These are interpreted to be remnants of the interaction 
withM81 approximately lOOMyrs ago |(|Yun et al.|1993a|l| 

While these kinematic parameters are confirmed by a large 
number of observations and studies, the contribution of the mag- 
netic field to the outflow is still a matter of debate. For normal 
spiral galaxies the mean field acu-dynamo in combination with 



galactic winds driven by supernova explosions in the spiral arms 
is the favoured explanation for the observed magnetic field pat- 
terns ( Krause 2009 , |Beck et al.|19 96 ) However, the role of the 
magnetic field in starburst galaxies is completely unclear. It can 
be a passive one, where the magnetic field is simply dragged 
along the velocity gradient of the wind, or an active one, where 
the magnetic field influences the gas motions in the outflow and 
it becomes a wind driven by cosmic rays and its associated field. 
The total magnetic field energy in normal spiral galaxies con- 
tributes only negligibly to the energy budget, but the high parti- 
cle densities and possible high magnetic field strengths in star- 
bursts can increase it to a value where it needs to be accounted 
for. A detailed study of these parameters is needed to constrain 
the inputs for simulations to understand material circulation and 
the observed morphology of galactic halos. This has a signifi- 
cant influence on models for the origin of magnetic fields and the 
magnetisation of the intergalactic medium in the early Universe, 
where starbursts and interactions are known to be more common. 
Starbursting dwarf galaxies in the early Universe can explain the 
magnetisation and enrichment of the intergalactic medium with- 
out a primordial seed field (Kron berg et al.|1999|) i 

Current instruments lack the sensitivity as well as the resolu- 
tion for a study of such objects in the early Universe, so that an 
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extrapolation from nearby observable targets is needed. Due to 
its proximity and large apparent size of Ilf2x4'3 and the avail- 
able data at nearly all wavelengths, M82 is the ideal target for 
such an analysis. In addition, its interaction scenario makes it 
possible to examine the influence of the group medium on a 
galactic superwind scenario hosting a starburst. 

|Yun et al.|fT994] l|and piynoweth et aT1 ( [2008| found a large 
amount of H i filling the whole group as well as three distinct 
spurs, which connect the main galaxies of the group: M81, M82, 
and NGC3077. Additionally, M82 is well know n for its large 
fraction of molecular gas in the central region (Aladro et al.l 
20TT]|Naylor et al.|2010]> Observations by|Greaves etaLlpOOO) 



Table 1. Basic properties of M82 



show a giant magnetic molecular bubble possibly blown out by 
the superwind. But how far this gas distribution influences the 
evolution of the whole group and whether it is reaccreted to 
their host galaxies is still not clear. An extensive star formation 
in the recently formed spurs could have enriched the intragroup 
medium with metals and magnetic fields ( |de Mello et al.|2 008 ) 
In how far the superwind in M82 is constantly fuelling the 
halo with fresh material is a matter of discussion. Evolution 
scenarios reach from an increased star-formation period lasting 
more than lOMyrs to the occurrence of several episodes of in- 
creased star-formation over the last Myrs. A detailed NIR-MIR 
spectroscopy by Forster Schreiber et al. (2003) suggests two suc- 
cessive starburst periods about 10 and 5 Myrs ago. This would 
explain the enrichment of the halo with ionised polycyclic aro- 
matic hydrocarbons (PAH) observed by Kaneda et al. ( 2010| )| 
which would otherwise have been destroyed by the starbursts vi- 
olent radiation. Additionally extended emission has been found 
~llkpc to the north of the galaxy in Ha, X-rays, and UV 
i Pevine & Bally|1999|phnert et al.|1999||Hoopes etaLp 005) 
which is spatially not connected to the main halo of M82. This 
might be a collision of the hot superwind with a preexisting 
neutral cloud (Lehnert et al. 1999) and a remnant of a previ- 
ously stronger star-bursting activity blowing material even fur- 
ther into the halo. The detection of hard y-ray emission by the 
| VERITAS Collaboration et al.| ( [2009] > further verifies a closed 
environment scenario since it favours the idea of M82 being 
a proton calorimeter ( Lacki et al.||201l| Additionally the high 
particle and energy densities in the core region lead to a sce- 
nario where even charged particles accelerated by supernova ex- 
plosions would not be able to escape this region and could not 
build up a synchrotron halo assuming the recent activity of star- 
formation and particle densities. 



Previous studies by Seaquist & Odegard ( 199Tj» and |Reuterj 
|et al.| ([1992} found emission, which was mostly confined to the 
core region of the galaxy and discussed the possibilities for reac- 
celeration and/or dynamo processes ( |Reuter et al.||1994] l inside 
the turbulent medium of M82, but could not sufficiently explain 
the diffuse emission reaching up into the halo. The improved 
data reduction and instrumental techniques together with more 
precise measurements over the whole wavelength regime now 
allows to study this prototypical starburst galaxy with increased 
detail. 

Since the synchrotron emission from gyrating electrons is the 
best tracer for magnetic fields, a complementary study at several 
different radio wavelengths between A3 cm and A92 cm is needed 
to get a full picture of the several loss processes and the mag- 
netic field morphology of this galaxy and its surrounding halo. 
The high flux difference between the highly luminous core and 
the very diffuse halo structure especially at longer wavelengths 
makes such a study challenging for the dynamic range require- 
ments of today's available instruments and data reduction tech- 
niques. Therefore data from the VLA at A3 cm and A6 cm and 



Name 

Alternative Names 

RA 2 ooo 

DEC2000 

Type 

Apparent Size (D 2 s) 
Distance 
Parallactic angle 
Inclination 



M82 

NGC3034, 3C231, Arp337 
09^55^52^7 
+69°40'46" 
10; Starburst 11 
11:2 x4:3 d 
3.52 Mpc b 
68 oc 

79 oc 



Tabl e 2. a NED ^Jacobs et al.| ( |200 9) c HyperLeda jHuchtmeier 
|& Skillman| ( [T998l l 



data from the WSRT at A22 cm and A92 cm were used for the 
following study. Both telescopes are well known for their excep- 
tionally good properties at these mentioned wavelengths. 

The goal of this paper is to investigate the outflow mecha- 
nism and to describe a revised scenario for the built-up of the 
observed halo. As a first step we investigate the magnetic field 
strength in M82 using an energy equipartition ansatz. These re- 
sults can then be used to calculate the different loss timescales 
and designate the dominant physical process influencing the 
wind in the two different regions. Arguments for the existing 
halo will be discussed and a scenario with a low filling factor of 
the ionised well as magnetised medium in the core region is de- 
scribed. These results and the comparison with observations at 
other wavelengths will be used to understand the outflow mech- 
anism and describe a revised scenario for the build-up of the 
observed halo. 

This paper is organised as follows: Section [2] shows the re- 
duced data and explains the special data reduction techniques 
used to reach the needed high dynamic ranges for this study. 
In section[3]the total intensity maps and vertical scaleheights are 
presented, section[4]shows the spectral index maps. Estimates for 
the magnetic field strength are given in section [5] and a detailed 
discussion about the different cosmic ray losses is presented in 
section [6] We summarise our outcomes in section [7] and discuss 
them in [8] 

2. Observations and data reduction 

2. 1. VLA observations 

The A3 cm and A6 cm data were obtained from the archive (al- 



ready published by Seaquist & Odegard ( 1991 ) and Reuter et al. 
(1992)), combined and the calibration improved for our stud- 
ies. We chose only the data from the compact D-configuration 
to match the resolution to our WSRT data and to restore most of 
the diffuse halo structure of M82. 

All observations used two subbands each at central frequen- 
cies of 8.465 GHz and 8.415 GHz and 4.885 GHz and 4.835 GHz 
respectively. Each subband had a bandwidth of 50 MHz, re- 
sulting in a total bandwidth of 100 MHz for each of the ob- 
served wavelengths. The observations were accompanied by cal- 
ibrator observations of 3C138, 3C286, 3C48, 1044+719, and 
0836+710 to determine the fluxes, phase corrections, polari- 
sation leakages and polarisation angles over the time of the 
observations. Data reduction and flagging was done using the 
Common Astronomy Software Applications package (CASA). 
Initially calibrator gains were applied independently to each sub- 
band and observation using the flux density scale given by Perley 
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Table 3. Observation parameters 



Wavelength 
(cm) 


Telescope 
and con- 
figuration 


Total 
bandwidth 
(MHz) 


Channels 

per 
subband 


Channel 
bandwidth 
(MHz) 


Subband 
frequencies 
(MHz) 


Polarisation 
products 


Total inte- 
gration time 
(min) 


Observing 
date 


3 


VLA-D 


100 


1 


50 


8465, 8415 


RR, LL, 
RL,LR 


160 


1990/01/04 


6 


VLA-D 


100 


1 


50 


4885, 4835 


RR, LL, 
RL,LR 


160 


1990/01/04 


6 


VLA-D 


100 


1 


50 


4885, 4835 
1753, 1737, 


RR, LL, 
RL,LR 


225 


1986/01/26 


1 o 


W oKl 


1 AH 


A/1 


n ^ i ^ 

\J.D 1 J 


1721, 1705, 
1689, 1673, 
1657, 1641 

1422, 1406, 


XX, YY, 
XY, YX 


3 AO 


ZUUj/UO/ 10 


11 
ZZ 


WCDT 
W oKl 


1 AO 
1 DU 


04 




1390, 1374, 
1358, 1342, 
1326, 1310 

1450, 1424, 


XX, YY, 
XY, YX 


7 AH 


ZUUj/Uj/ 10 


22 


WSRT 


80 


128 


0.078 


1419, 1398, 
1374, 1349, 
1326, 1303 

321, 325, 


XX, YY 


720 


2003/02/09 


92 


WSRT 


20 


128 


0.019 


332, 336, 
372, 376, 
380, 385 


XX, YY 


720 


2003/04/30 



and Butler (2010). Then self-calibration was performed in the 
same way calibrating four times on phase only followed by one 
simultaneous amplitude and phase calibration. For each itera- 
tion a model was derived using the Clark-Clean algorithm (Clark 
1980) Inside each self-calibration loop the clean iterations were 
increased to include more of the diffuse flux in the model and the 
time interval for solutions decreased until the integration time 
was reached. The data were then cleaned down to a 3<x-noise 
level and baseline-based corrections were applied using a long 
solution interval of 10 minutes for the total power images to in- 
crease the dynamic range. Integrated fluxes for the total power 
images before and after the baseline-based corrections were ap- 
plied differing by less than 1.5 %. Not enough flux was available 
in the Stokes Q and U images, so that the baseline-based cor- 
rections were not used for the production of the polarisation im- 
ages. For each clean step masks were used from previous data re- 
ductions iteratively including more sources inside the field. The 
resulting final images of each observed wavelength were com- 
bined using an inverse square weighting followed by primary 
beam correction. 

This subband-based calibration strategy lowers the fre- 
quency dependence of the gain solutions as well as the flux vari- 
ations of sources inside the field over the complete bandwidth. 
Care has to be taken for observations with low flux intensities 
since the described method lowers the signal-to-noise for the de- 
termination of the gain solutions, which can lead to an image 
quality worse than calibrating using the whole data in all. 



2.2. Westerbork observations 

The WSRT archive was inspected for archival datasets of M82. 
One dataset at 18 cm, two at 22 cm and one at 92 cm were found. 



The 18 cm data and one of the 22 cm datasets were recorded 
in all four cross-correlations. Therefore these data contained 
polarisation information. For the other only the XX- and YY- 
correlations were available. From this data only a total power 
image of the 18 cm dataset has been published so far by |Braun 
et al.|l|2007)| 

All datasets were first read into the Astronomical Image 
Processing System (AIPS) to apply the system temperatures. 
Since AIPS can normally only handle circular polarised feeds, 
the data headers were changed to circular polarisation during 
this step and changed back to linear polarisation later. The data 
were converted into the CASA MS-format and inspected us- 
ing RFI GUI, the graphical frontend to RFI console (Offringa 
et al. 2010) A special flagging strategy was developed for each 



wavelength and source (calibrators/M82). Since the bandpass for 
the WSRT is very steep at the subband edges, a preliminary 
bandpass was applied to the data, which was not used for cal- 
ibration during later stages of the data reduction. For this the 
calibrator observation least affected by RFI (Radio Frequency 
Interference) was used for the appropriate wavelength. This pro- 
cedure helps RFI console to distinguish more easily between real 
RFI and bandpass effects. After this initial flagging each dataset 
was again checked for leftover RFI, which was then flagged 
manually. 

Calibrator gains were determined from the unpolarised cal- 
ibrator, and polarisation leakage and angle terms were deter- 
mined for the polarised calibrator. This was done for each chan- 
nel separately. The model parameters for each calibrator were 
based on a fourth order polynomial with coefficients A, B, C, and 
D from the VLA Calibrator ManuaQ Additionally the models 



http : //www . via . nrao . edu/astro/calib/manual/baar s . 



html 



3 



B. Adebahr et al.: M82 - A radio continuum and polarisation study 
Table 4. Summary of calibrator parameters 



Calibrator 


A 


B 


C 


D 


RM 

(rad/m 2 ) 


V 


$0 

(rad) 


3C138 (22cm) 
3C138 (18cm) 


1.0076 


-0.5562 


-0.1113 


-0.0146 


-1.2867 


0.0750 
0.0845 


-0.1379 


CTD93 


0.7110 


-0.0230 


-1.3190 


0.4938 











3C147 


1.4486 


-0.6725 


-0.2112 


0.0408 











3C295 


1.4674 


-0.7735 


-0.2591 


0.0075 












of the polarised calibrators were given a fractional polarisation 
p, an intrinsic polarisation angle <t>o and a rotation measure RM. 
This information was gathered from the newest EVLA measure- 
ments by R. Perley (priv. comm.). According to the equations 



j _ 2QA+B logv+C (logv) 2 +D (logv) 3 

Q= p- 1 ■ cos 



U 
V 



p ■ I • sin 
0. 



•((Do 



+ RM- 



RM 



1 

v) 



(1) 
(2) 

(3) 
(4) 



where I,Q,U, and V are the flux densities of the Stokes param- 
eters, c is the speed of light, and v the channel's frequency, and 
using the relations for the linear feeds of the WSRT 

XX = I + Q XY=U + V YX=U-V YY = I - Q (5) 

the fluxes for all four cross-correlations could easily be calcu- 
lated for each individual channel. A summary of the parameters 
used is listed in Tab. [4] 

This channel based derivation of the calibrator gains min- 
imises the effects of the 17MHz-ripple in the pointing centre 
of the beam and reduces the instrumental polarisation signifi- 
cantly, which can vary by up to 1.5% across the frequency band 
(Brentjens 2008) Fig[T] shows the difference between assuming 
a constant gain across one subband and the channel-based cali- 
bration performed here. 

After this initial calibration step, each dataset was self- 
calibrated with different strategies described in the next subsec- 
tions. 



2.2.1 . 1 8 cm and 22 cm datasets 

The first 22 cm dataset contained a single 12hr run with a to- 
tal bandwidth of 80 MHz split into eight subbands centered at 
1450, 1424, 1419, 1398, 1374, 1349, 1326, and 1303 MHz. Each 
subband was covered by 128 channels in two cross-correlations. 
The integration time was 30 seconds, which was sufficient to 
sample the ionospheric fluctuations. The observed calibrators for 
this dataset are 3C147 and 3C295. The second dataset contained 
a complete 12hr-run with frequency switching between 18 cm 
and 22 cm every 5 minutes yielding an integration time of ~ 360 
minutes for each frequency band. The subbands were centred 
at 1753, 1737, 1721, 1705, 1689, 1673, 1657, and 1641MHz 
for the 18 cm-band and at 1422, 1406, 1390, 1374, 1358, 1342, 
1 326, and 1310 MHz for the 22 cm-band. Each subband was cov- 
ered by 64 channels in all four cross-correlations with an inte- 
gration time of one minute. The observed calibrators for these 
two datasets are 3C138 and CTD93. 

After applying the calibrator gains with the strategy de- 
scribed in Sect. |2.2| both datasets were exported to MIRIAD 



8.8 
8.7 
8.6 
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8.4 
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Fig. 1. Flux of the central pixel of the used calibrator 3C138. 
The top plots show the flux in Stokes I and the bottom ones the 
flux in Stokes Q, U, and V. The improvement between the sub- 
band based calibration (left) and the channel based one (right) is 
obvious. 



for faster imaging and cleaning. Due to the extremely luminous 
extended core of M82, a special self-calibration technique, de- 
scribed in Sect. 2.2.2 had to be used to reach a dynamic range 
of ~50000. For each step in this data reduction involving clean- 
ing masks were used from previous data reductions to include 
all visible flux inside the field. 

Because of the better data quality of the observations with 
80 MHz bandwidth, the shorter integration time and the simi- 
lar subband frequencies, clean models from the former dataset 
at appropriate frequencies were used to start the self-calibration 
process for the dataset with 160 MHz bandwidth at 22 cm. This 
leads to 13 usable images, which were then combined with an 
inverse square weighting and primary beam corrected. 

The 18 cm dataset was only used for polarisation purposes 
and therefore no total power images were produced. 



2.2.2. Calibration technique for strong extended sources 

The East-West alignment of the WSRT and the observation 
mode with only one calibrator observation in the beginning 
and one at the end makes self-calibration the most important 
step in the data reduction for high dynamic -range images. One 
has to make sure that especially the first models during the 
self-calibration process do not include spurious features, which 
would then be transported into images and models for later 
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stages. These features are normally easily visible in the images 
since they are biased towards calibrator observations and array 
layout. This means for the strong extended source M82 that 
artificial clean components would be used which elongate the 
source in the East- West direction. 

Unresolved sources are normally very good calibrators since 
their total emission is covered by a clean beam. This makes 
them much better sources for self-calibration than strong ex- 
tended sources since their clean component model is only one 
delta-function and not multiple delta-functions resembling the 
complex diffuse emission like it is for M82-type sources. The 
clean components for these sources would resemble the artifi- 
cial elongation in the East- West-direction. Therefore, the best 
choice would be to construct a (u,v)-dataset of the M82-field 
with a dominating flux contribution from point sources. This was 
achieved by using the following technique: 



1. The channel-based initial calibration (Sect. 2.2 1 was applied 



to the dataset. This minimises the contribution of the 17 MHz 
ripple in the central part of the field. 

2. We applied a normal self-calibration approach for the entire 
field using only baselines > 4 kA for the first self-calibration 
loop with a high solution interval. This parameter was re- 
duced as well as the solution interval during further loops. 

3. We subtracted M82 from the dataset using the gains and the 
clean model from the last self-calibration loop. Afterwards 
the inverted gains from the same loop were applied resulting 
in a dataset without M82 and only the gains from the initial 
calibration applied. 

4. A self-calibration on this dataset was applied using only the 
point sources yielding an average gain over the entire field. 
After two rounds of self-calibration with long solution inter- 
vals and (u,v)-ranges > 4kA the gains were transferred to 
the original initially calibrated dataset containing all point- 
sources and M82. Then the self-calibration was applied to 
this dataset in the normal manner reducing the solution inter- 
val and including smaller (u,v)-ranges to improve the model. 

This technique enables the use of the high-signal-to-noise- 
ratio of the strong extended source without artificially including 
clean components from the East- West-biasing. A comparison of 
the image quality before and after the described technique was 
applied is shown in Fig. [2] 

2.2.3. Polarisation calibration 



Due to high rotation measures already detected by |Reuter et al.| 
( |1994| the derivation of a model for all four Stokes parame- 
ters for each single subband is not sufficient. The Stokes Q and 
U fluxes vary significantly over the bandwidth of one subband. 
This would result in a wrong model for a significant number of 
channels in each subband, shifting flux from Q to U or vice versa 
and producing artificially false rotation measure components. 
Therefore, self-calibration for each single channel was done in- 
dependently as described by de Bruyn & Brentjens (2005) The 
high flux of M82 was sufficient to reach the signal-to-noise re- 
quired for a satisfactory dynamic range in the Stokes Q- and 
U-parameters. 

The whole calibration was done in MIRIAD using a script 
and including all (u,v)-ranges in the self-calibration. Masks from 
the total power data reduction were used to include all diffuse 
emission in the model. The final self-calibration loop used an 
amplitude and phase calibration with a solution interval of one 
minute. The final images were then cleaned to two times the 



noise level and primary beam corrected. After that all images 
were inspected for significant artefacts due to calibration errors 
and leftover RFI by eye. At the time of the observation the sun 
was at a maximum of its activity, which led to significant signals 
in the polarisation images far off the pointing centre. Since these 
signals are heavily time- and frequency-dependent, a subtraction 
of the signal was not possible, so that these images had to be 
discarded, but images at frequencies with a lower contribution 
of the signal from the sun were usable. Additionally any images 
suffering from a significant elongation of the beam as a result 
of time-based flagging of the data of an east-west-interferometer 
like the WSRT and images with a noise floor three times higher 
than the average were not used. This yielded only 396 (187 for 
the A18cm and 209 for the A22cm observations) out of 1024 
final images for each polarisation product. 

2.2.4. 92 cm dataset 

The 92 cm dataset contained a single 12 hr exposure with a total 
bandwidth of 20 MHz. Each of the 8 subbands was sampled by 
128 channels in two cross correlations. The central subband fre- 
quencies were set to 321, 325, 332, 336, 372, 376, 380, and 385 
MHz. The observed calibrators were 3C147 and 3C295. The in- 
tegration time was 30 seconds, which is sufficient to sample the 
ionospheric fluctuations during the night. Unfortunately some 
part of the observation was carried out during the day and at 
sunset, where the interference from the sun is especially high on 
the short baselines. Subtracting a clean model of the sun after the 
initial calibration was not successful, so that several hours on the 
shorter baselines had to be flagged. 

The data was then exported to MIRIAD and self-calibration 
was carried out for each subband independently using masks 
during the clean step to include all diffuse emission in the field. 
The small bandwidth of each single subband of only 2.5 MHz 
made it possible to use a normal Clark-clean algorithm rather 
than a multi-frequency deconvolution (Sault & Wieringa 1994) 
to reach a dynamic range of ~ 10000. The subband images were 
then combined with an inverse square weighting and primary 
beam corrected. 



3. Total power emission 

3.1. Morphology 

The maps showing the total power emission are presented in Fig. 
[3] Fig . [4j Fig . [5] and Fig . [6] at wavelengths of 3 cm, 6 cm, 22 cm, 
and 92 cm respectively. 

The morphology of the total power emission depends 
strongly on the observed wavelength. While the maps at 3 cm 
and 6 cm show only emission in the central disk and some extra- 
planar emission to the north and south of the very central region 
of M82, the maps at 22 cm and 92 cm show the structure of the 
radio halo and the disk emission. Overall the halo emission in 
the southern part seems to be less extended than in the northern 
part, but the whole halo emis sion correlates wi th the Ha- as well 
as the PAH + -morphology ( |Kaneda et al.|2010[) | 

The most obvious difference between the northern and south- 
ern halos are two extensions accompanied by a suppression of 
emission in between them in the northern halo. While the 92 cm 
map shows basically the same features towards the north and 
south with a smaller extent towards the south, the 22 cm mor- 
phology is completely different in the southern halo. Here the 
morphology is dominated by an extension south of the starburst- 
ing core with lacking emission in the southwestern part. 
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Fig. 2. A comparison of the image quality using a subband-based initial calibration technique and standard self-calibration procedure 
(left), using a subband-based initial calibration and the self-calibration technique described in Sect. 2.2.2 (middle), and using a 



channel-based initial calibration technique and the self-calibration technique described in Sect. 2.2.2 (right). The biasing of the flux 
towards the East- and West-direction and resulting sidelobes are easily visible in the first image. These effects nearly vanish in the 
second image and completely disappear in the last one. The images show only the first subband of the observations with 80 Mhz 
bandwidth at A22 cm without any primary beam correction. 




Fig. 3. Total power radio continuum contours at A3 cm from 
the VLA-observations overlaid on a Ha image from the WIYN 
3.5 m Observatory provided by M. Westmoquette et al. (priv. 
comm.). Contours start at a 3<x level of 50 |aJy/beam and increase 
in powers of 2. The beam size is 7'.'6x7'.'3 and is shown in the 
bottom left corner of the image. 



Fig. 4. Total power radio continuum contours at A6cm from 
the VLA-observations overlaid on a Ha image from the WIYN 
3.5 m Observatory provided by M. Westmoquette et al. (priv. 
comm.). Contours start at a 3<x level of 120 nJy/beam and in- 
crease in powers of 2. The beam size is 12'.'5xl 177 and is shown 
in the bottom left corner of the image. 



While the extent of the southern extensions at A92cm is 
about the same, it is different in the northern halo at A22cm 
as well as at A92 cm. Here the northwestern spur seems to be 
more extended than the northeastern. McDonald et al. (2002} al- 
ready mentioned a difference in the distribution of compact radio 
sources in the core region between the eastern and western part. 
The western part seems to be denser populated by supernova 
remnants than by HII regions while it is vice versa in the east- 
ern part. This could mean that the western part is further evolved 
and ejected the material into the halo earlier than the eastern part 



where a higher supernova activity is expected within the next 
Myrs. 

Remarkable is the morphology of the radio emission along 
the disk of M82. Obser vations by |Seaquist & Odegard| (|1991|l 



and Reuter et al. 



( 1992| already showed this extent of the emis 



sion. Therefore it is most likely that the radio morphology of 
M82 consists of a superposition of a disk and a biconical out- 
flow. 

While the lowest contour levels for the 3 cm, 6 cm, and 22 cm 
maps have similar noise levels, the large halo emission is only 
visible in the 22 cm and 92 cm map. This might be partly due 
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Fig. 5. Total power radio continuum contours at A22 cm from the 
WSRT-observations overlaid on a Hoc image from the WIYN 
3.5 m Observatory provided by M. Westmoquette et al. (priv. 
comm.). Contours start at a 3<x level of 90 |jjy/beam and increase 
in powers of 2. The beam size is 12'.'7xll'.'8 and is shown in the 
bottom left corner of the image. 
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Fig. 6. Total power radio continuum contours at A92 cm from the 
WSRT-observations overlaid on a Ha image from the WIYN 
3.5 m Observatory provided by M. Westmoquette et al. (priv. 
comm.). Contours start at a 3cr level of 1.8mJy/beam and in- 
crease in powers of 2. The beam size is 43'.' 1 x3976 and is shown 
in the bottom left comer of the image. 

Table 5. Integrated fluxes for our M82 observations 



A (cm) 


Flux (Jy) 


3 


2.20 ±0.11 


6 


3.23 ±0.16 


22 


7.41 ±0.37 


92 


13.83 ±0.69 



to lower sensitivities of the A3 cm and A6 cm maps to extended 
emission because of missing spacings. Comparing our integrated 
flux of 3.23 Jy + 0.16 Jy with the single dish observations of 
Gregory & Condon (1991) at A6cm, who measured a flux of 
3.96 Jy + 0.60 Jy, may support this statement. The surrounding 
point-sources were not subtracted for the single-dish measure- 
ment due to the lacking resolution, but adding up the flux of 
all point-sources in our interferometer image within the beam 
of the single dish telescope yields only ~ 20mJy. The largest 
observable scale of the VLA at A3 cm is 5'A and 9'at A6 cm in 
comparison to the WSRT with 21' at A22 cm and 81' at A92 cm. 
This is always significantly larger than the observed structures at 
the appropriate wavelengths. Additionally, the VLA in its used 
compact D-array configuration has a lot of short baselines, which 
make it very sensitive to structures larger than the ones observed 
with an extent between 1!5 ~ 2'. A comparison of our obser- 
vations with the ones from Soid a et aL] ( |201 1[ ), who also used 
the VLA D-array configuration for their studies of NGC5775, 
is more informative. They could recover emission with a size of 
3'.5 at A3 cm and A6 cm. A combination of their VLA measure- 
ments with single-dish data from the Effelsberg telescope only 
recovered 3 % more flux at A3 cm and no increase was recog- 
nised at A6 cm. The rapid decrease of the emission in direction 
away from the galaxy is stronger than expected from missing 
spacings alone. Therefore we conclude that we can neglect in- 
fluences of missing spacings to our flux values. 

3.2. Total intensity distribution 

To isolate the emission of the galaxy from contributions of back- 
ground sources, all point-like sources exceeding three times the 
rms noise in the images were fitted with Gaussians in the im- 
age plane and subtracted. This method was successful for all 
sources except the double source to the southeast of M82 due 
to its complicated structure and the fact that it is embedded in 
the diffuse emission. Therefore this source was left in the im- 
ages and for parts of the data analysis where its flux contributes, 
the data points were not included in the analysis. 

3.2.1. Radial total intensity profile 

The resulting total intensity profiles along the major axis are 
shown in Fig. [7] for all four observed wavelengths. The indi- 
vidual values were obtained by regridding the total power maps 
to a SLOAN G' image of the galaxy and integrating the flux 
over each pixel column. Since we are only interested in the over- 
all shape of the profile and the relative differences between the 
wavelengths, the absolute integrated fluxes are irrelevant. 

While the emission in the central part is dominated by the 
star-forming nucleus of M82, the 22 cm- and 92 cm-observations 
show significant emission in the disk of the galaxy. Remarkable 
is a slight shift of the emission maximum to the west visible in 
the higher frequency images due to their better resolution. This 
seems to be a consequence of the higher star-forming activity 
in the western part of the central region of the galaxy (Muxlow 
et al. 1994) and might even cause the difference between the ex- 
tent of the spurs in the northern region. The most extended fea- 
ture reaching up to ~4.5 kpc into the halo is positioned directly 
north of this region. Here the cosmic ray gas may be able to es- 
cape the galaxy easier since the increased star-formation simply 
blew away the dense medium providing a channel towards which 
the gas can escape from the central starburst. 
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with the complementary error function: 
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Fig. 7. Integrated total intensity for each individual pixel column 
plotted versus the distance from the centre of the galaxy along 
the major axis. The flux was integrated over the total power 
emission along the minor axis of the galaxy. The maximum at 
-100" is the consequence of the source left inside the diffuse 
emission (see Sect. 13. 2b. 



3.3. Vertical scaleheights 

The observed radio emission of an inclined galaxy is a mixture 
of the halo emission and the projected disk emission contribut- 
ing to the halo emission (Du mke et al.|1995] > In order to simu- 
late the latter, we compressed the intensity distribution along the 
major axis by the inclination of the galaxy. The resulting distri- 
bution was then convolved with the HPBW of the observation to 
account for beam smearing effects. The FWHM of a Gaussian 
fitted to the major axis profile of this image is then the effective 
beam size. This led to values of 16"(265 pc) at A3 cm, 19" (330 
pc) at A6 cm, 2 1 " (350 pc) at A22 cm, and 50" (860 pc) at A92 cm. 
All emission outside this distance from the major axis can be at- 
tributed to extraplanar emission. 

To analyse the propagation of cosmic rays and their losses 
a parameter, which is independent of the total flux, is needed. 
Therefore the vertical scaleheights were used. This approach 



has already been performed for other objects (Dumke & Krause 
[19981 |Heesen et al.|2009||Soida et al.|20TT] l and makes the val- 
ues comparable to each other. 

According to|Dumke et al. ( 1995 i the fitted function in case 



of a Gaussian distribution is 



W,„ 



s (z) = 



WQZ0 



: exp 



2ct 2 +z 2 q ) 



(6) 



with cr being the effective beam size, wq the maximum of the 
distribution and zo the scaleheight. In case of an exponential dis- 
tribution this becomes 

_2 



W exp (z) 



wo 



exp 



z 
"2^2 



/cr 2 - ZZ0 . . 

exp — — | erf c 



\ V2, 



o-zo 



+ exp 



cr 2 + ZZQ V 



V2crz ) erf ° I V2crzo ) 



cr - ZZQ 

V2o-z 

cr 2 +ZZo\ 



(7) 



erfc x = 1 — erf x ■ 



2 



exp(-r )dr 



(8) 



Integrations over rectangular areas with a width of 60" parallel 
to the galaxy's minor axis were used as an input to a Gaussian, 
a one component exponential and a two component exponential 
fit. This divides the galaxy into 3 strips each north and south of 
the major axis for the A22 cm and A92 cm data and one strip each 
for the A3 cm and A6 cm data resembling a good database for the 
spurs in the east and west and the outflow cone in the centre. 
This integration technique ensures that the fits do not suffer from 
small scale fluctuations in the images. 



We used the fitting routines by Dumke et al. ( 1995) which 



take the calculated effective beam size into account. Our mea- 
sured values and the three different types of fit functions are 
presented in Figures [8] to 1 1 the deprojected scaleheights are 
presented in Table [6] 

Data points where the background radio source in the south- 
west contributed to the diffuse emission were not used for the 
fitting routines. Errors cr s for the data points were calculated us- 
ing: 

= J^Ls + °\> (9) 



where cr rms is the error due to the noise in the image, which 
was estimated by selecting an emission free area close to the 
centre of the image and calculating the standard deviation, and 
°"b = Vms V" me calculated baselevel error in the image, with 
n being the number of integrated beam areas. Additionally the 
standard deviation inside each integration area was calculated 
and in cases where it was higher than the calculated baselevel 
error, <x B was replaced by this value. This is especially critical in 
areas with strong gradients like the core region, where the actual 
fluctuations of the source over the area is dominating the error 
for the integration. 

3.3.1. Scaleheights at A3cm and A6cm 

Since the observations at A3 cm and A6 cm only show emission 
in the central part of the galaxy, the analysis of the scaleheights 
is mostly limited to the disk and inner halo emission. Fig.[8]and 
Fig.|9]clearly show that the one component and two component 
exponential fits fit much better than the Gaussian fit. 

The absolute values for the A3 cm data and A6 cm data do 
not differ a lot from each other. This is a clear sign for a homo- 
geneous halo property at these wavelengths for this part of the 
galaxy. Therefore an average value of the scaleheights can be 
used by simply calculating the mean. This yields average scale- 
heights of 20 ± 5 pc for the thin disk and 140+30 pc for the thick 
disk at A3 cm and A6 cm. 

3.3.2. Scaleheights at A22cm and A92cm 

Looking at the emission profiles at longer wavelengths with high 
sensitivity shows more extended disk and halo emission. The 
fits in Fig. [TUJand Fig. [TT] clearly show that there is no function 
which fits all the data points, but for most cases the two com- 
ponent exponential fit is the best solution. The exception is the 
southwestern part of the galaxy, where the emission drops very 
rapidly in the z-direction. This seems to be a clear sign of high 
losses of the cosmic ray electrons in this area. No successful fit 
could be found for the northwestern part at A22 cm. 
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Fig. 8. Intensity distribution and its errors for the A3 cm-data. 
Each data point reflects an integration over an area of 60"x4.5". 
The blue line shows a Gaussian fit, the green line a one compo- 
nent exponential fit, and the red line a two component exponen- 
tial fit. The parameters for the fits are listed in Tab. [6] 
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Fig. 9. Intensity distribution and its errors for the A6 cm-data. 
Each data point reflects an integration over an area of 60"x7". 
The blue line shows a Gaussian fit, the green line a one compo- 
nent exponential fit, and the red line a two component exponen- 
tial fit. The parameters for the fits are listed in Tab. [6] 

Remarkable is the difference in the scaleheights between the 
northern and southern part of the galaxy. Using again the mean 
of the A22cm yields 690 ± 50 pc for the northern thick disk 
and 390 ± 70 pc for the southern one, and for the A92 cm emis- 
sion 820 ± 170pc and 450 ± 20pc respectively. Since most 
of the galaxies in the M81-M82-group core are located in the 
south of M82, the group medium is expected to be denser to- 
wards this direction. The reduced scaleheights in the southern 
part are a good argument for higher losses of cosmic rays to- 
wards the group centre. [Kaneda et aLj pOlO) detected a stronger 
radiation field only in the southeastern part of the galaxy which 
favours higher inverse Compton losses in this area, but this can- 
not explain the overall reduced extent of the southern halo and 
its reduced scaleheights. A slower southern outflow would be a 
possible explanation, but this would need a lower starburst ac- 
tivity and/or density in the southern central part to explain the 
reduced kinetic energy input. Arguing for a motion of M82 to- 
wards the group centre and a compression of the cosmic ray gas 
as has been seen in clusters of galaxies cannot be confirmed by 
simulations ( |Yun et al.| 1993a) , but is also not fully ruled out due 
to the uncertainties in the possible evolution scenarios. 

Since it is most likely that secondary electrons are produced 
in the core of M82 and its surrounding region up to several hun- 
dred kiloparsecs (see Sect. [6] for a detailed discussion) a one 
component exponential fit was also used to model loss processes 
that are independent of the core region. No strong difference is 
seen between the scaleheights at A22 cm and A92 cm. Averaging 
the scaleheights for the thick disk of both wavelengths results in 
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Fig. 10. Intensity distribution and its errors for the A22 cm-data. 
Each data point reflects an integration over an area of 60"xl5". 
The blue line shows a Gaussian fit, the green line a one compo- 
nent exponential fit, the red line a two component exponential 
fit, and the cyan line a one component exponential fit to the halo 
only. The parameters for the fits are listed in Tab. [6] 



Table 6. Scaleheights between A3 cm and A92 cm for all success- 
ful fits. The indices in the first column show the position of the 
fit. Zo,2e,c and Zi,2e,c are the scaleheights for the first and second 
component to the two-component exponential fit and Zo,u,h i s the 
scaleheight of the one-component exponential fit to the halo. 



Strip 


A(cm) 


Zo>, r (pc) 


Zl,2e,c(P c ) 


zo,w,(pc) 


x s,c 


3 


16 


121 




X N,C 


3 


22 


181 




X S,C 


6 


22 


134 




X N,C 


6 


20 


129 




X S,E 


22 


72 


445 


819 


x n,e 


22 


146 


663 


874 


x s,c 


22 


40 


335 


441 


X N,C 


22 


143 


705 




xs.w 


22 


20 


334 


359 


X N,W 


22 








X S,E 


92 


114 


450 


811 


X N,E 


92 


107 


745 




x s,c 


92 


58 


471 


459 


X N,C 


92 


86 


700 




X S,W 


92 


431 


431 


541 


Xff,W 


92 


77 


1017 





610 + 230pc. All values for the thick disk in Tab. |6]are still by 
a factor of 2-4 smaller than the ones from normal spiral galaxies 
flDumke & Krause|1998||Krause|2011)| 
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Fig. 11. Intensity distribution and its errors for the A92 cm-data. 
Each data point reflects an integration over an area of 60"x30". 
The blue line shows a Gaussian fit, the green line a one compo- 
nent exponential fit, the red line a two component exponential 
fit, and the cyan line a one component exponential fit to the halo 
only. The parameters for the fits are listed in Tab. [6] 



4. Spectral index distribution 

4.1. Morphology 

Spectral index maps were produced between A3 cm and A6 cm 

and A22 cm and A92 cm 



(Fig. 
(Fig. 



12 1, A6 cm and A22 cm (Fig. [T3J 

14 1 using a 5o~-cutoff level for the total power images. To 
check the contribution of different (u,v)-coverages for the dif- 
ferent wavelengths, each spectral index map was once produced 
by limiting the used baselines to match the shortest and longest 
baselines of each other and once using all data and convolving 
the final maps with the same beam after image restoration and 
primary beam correction. No significant difference was visible 
so that the latter maps were used due to their higher sensitivity. 
Additionally error images were calculated leading to errors of 
cr a = 0.02 for the core regions and up cr a = 0.2 for regions close 
to the cutoffs. We use the convention / oc •/" for the following 
analysis. 

Thermal emission is known to contribute less than 10% of 
the total flux density of M82 at A6 cm and even does not dom- 



inate the spectrum at frequencies up to 25 GHz ( |Seaquist & 
Odegard 1991; Kl ein et al.| fl988), which implies a negligible 
thermal contribution to the measurements at A22 and A92 cm. 

All maps show an overall flat spectral index between a = 
-0.6 and a = -0.8 for the central part of the galaxy, which de- 
creases towards steeper values of a — - 1 .5 for the outer parts of 
the disk and halo emission. The average spectral index gradient 
in the halo is relatively flat with an index of a = —1.2 and does 
not steepen with increasing distance from the midplane (Fig. 
[14) . Two remarkable regions, one towards the northwest and and 
one to the south, are visible between wavelengths of A3 cm and 
A22cm, where the spectral index first steepens and then rises 
with increasing distance from the disk. Maps from Seaq uist &| 
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Fig. 12. Spectral index distribution between A3 cm and A6 cm. 
Contours start at -2.4 and end at -0.6 with an increment of 0.3. 
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Fig. 13. Spectral index distribution between A6 cm and A22 cm. 
Contours start at -2.1 and end at -0.3 with an increment of 0.3. 



Odegard ( 1991) who used different data for A22 cm and A92 cm, 
independently show the same characteristics at these locations. 

Looking into these two regions at spectral indices derived 
from longer wavelengths does not show this behaviour due to the 
beam size but following the same direction deeper into the halo 
shows spectral indices of a = -2.5. While the northern region 
is located directly in the outflow cone indicating high inverse 
Compton, adiabatic, and/or synchrotron losses in this region, the 
southern region can be associated with higher inverse Compton 
losses towards a recently formed OB-association due to the gen- 
erated stronger radiation field. 

4.2. Disentangling the core and halo emission 

The loss processes of the cosmic ray electrons are expected to be 
different in the high density region of the core and in the halo. 
To analyse this the flux at each wavelength was integrated once 
over the whole galaxy, once in the starburst region with a typical 
radius of 450 pc ( |Leeuw & Robson 2009 ) and once over the halo 
(Fig. p3). An additional measurement could be acquired from 
|Braun et~a l. (2007) for the core region, but the dynamic range 
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Fig. 14. Spectral index distribution between A22 cm and A92 cm. 
Contours start at -2.4 and end at +0.3 with an increment of 0.3. 



of their maps was not sufficient to extract a measurement for the 
halo. 

The spectral index is defined as 



S = S v° 



(10) 



The slope of the fits to the intensities S in Jy at the wavelengths v 
in GHz resembles the spectral index. Looking at the data points 
in a log-log-diagram (Fig. 15 i shows that only the halo compo- 
nent can be fitted with this function yielding »/, = -1.20 ± 0.01. 
Using only this function for the total emission resembles the av- 
erage of the spectral index over the whole galaxy with a a = 
-0.52 + 0.04, but this is heavily influenced by the core region, 
where Eq. 10 does not provide an adequate fit. 

A large amount of flux in the core region of M82 has been re- 
solved into ~ 100 independent discrete sources, which are a mix- 
ture of supernova remnants and ionised regions around newly 
forming star-clusters ( |Fenech et al.|2010||McDonald et al.f2 002 



Will s et al.| 1997| > The prior authors found indications that par- 



ticles in these dense regions are suffering from free-free absorp- 
tion and bremsstrahlung self-absorption, which reduces the ob- 
served flux for low frequencies. Under the assumption that the 
measured flux for the resolution of the presented observations is 
dominated by these objects in the core region, the opacity t of 
the medium can be included in Eq. 10 resulting in 



S = S v a e- T 
in case of ionisation losses and 

S =S ( )V 2 (l-0 
in case of bremsstrahlung self-absorption with 
8.2 ■ 10~V 21 EM 



(11) 



(12) 



T = 



-1.35 



where EM is the emission measure of the region in pc cm ° and 



T P the elec tron temperature of the warm medium in K (Wills 



et al. 



19971). Assuming T e = 10 4 K, the emission measure as 
well as the spectral index can be calculated by fitting Eq. [TT] 
and 12 to the data. Only the fit for the ionisation losses con- 



verged sufficiently, which suggests that these are the dominant 
losses at these wavelengths. This results in a spectral index of 
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Fig. 15. Integrated Flux of M82 for the whole galaxy, the core re- 
gion and the halo plotted vs. the frequency. The additional point 
at A18 cm is from Braun et al. (2007) The integrated emission 
over the whole galaxy shows a turnover in the spectral index, 
which can be attributed to free-free absorption in the core re- 
gion. The dashed lines represent the expected flux in case of no 
free-free absorption and the dotted line the best fit in case of 
bremsstrahlung self-absorption. 



a, = -0.67±0.01 for the total emission and of a c = -0.62±0.01 
in the core region with emission measures of EM, = 1.12 ± 
0.02 • 10 5 pccnT 6 and EM C = 3.16 + 0.1 • 10 5 pcctrT 6 respec- 
tively. Using the calculated values results in an opacity of the 
core region of t c = 0.92 at A92 cm. Since the emission measure 
is defined as 

Jl. 



EM 



K 





ds, 



(14) 



where n e is the electron density of the ionised medium and so 
the line-of-sight through the medium, n e can easily be calcu- 
lated. Using .so = 900 pc results in an electron density of n e = 
18.7 cm 3 . Assuming that nearly all electrons are confined inside 
the star-forming regions and supernova remnants, the filling fac- 
tor f e can be calculated. Using a density of n e = 1000 cm 3 from 
Westmoquette et al. (2009) for the star-forming regions leads to 
a filling factor of f e = 1 .9%, which is confirmed by theoretical 
predictions from Youst-Hull et al. (priv. comm.). 

These calculations suffer from a contribution to the inte- 
grated area of the disk emission in the front and on the back- 
side of the galaxy as well as the unknown morphology in the 
core region. But anyhow they are good enough for an order of 
magnitude estimate. 



(13) 5. Magnetic field strength 



Assuming equipartition between the magnetic field energy and 
the energy density of the cosmic rays, one way to determine the 
total magnetic field strength is using the revised equipartition 
formula by Beck & Krause (2005 ) 

Due to the complicated structure of M82, the high dynamic 
ranges inside the galaxy in the total intensity and the complex 
environment it is embedded in a pixel based method for the de- 
termination of the magnetic field strength and cosmic ray losses 
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was used. This ensures that the steep gradient over the galaxy of 
the total power emission is accounted for and the error for the 
calculation of the magnetic field strength due to the integration 
over the line-of-sight is minimised. The A22 cm fluxes were used 
for the entire calculation because of its superior dynamic range 
over the A92 cm map. 

A cylindrical geometry was assumed for the synchrotron 
emitting regions, for which the projected profile is a Gaussian 
with constant values along the minor axis, and a Gaussian profile 
along the major axis. The full width half maximum for these ge- 
ometries was calculated by convolving the 22 cm map with the 
92 cm beam, to match the resolution to the spectral index map 
making the result comparable, integrating the emission along the 
minor axis and fitting a Gaussian to the resulting profile yielding 
a FWHM of 930 pc. Moreover the assumed non-thermal spectral 
index was used from the halo fit in Fig. 15 since the equipartition 
in the formula is assumed to take several Myrs. This leads to reli- 
able magnetic field strengths only in the halo region. To estimate 
the magnetic field strength in the core region a spectral index of 
a = -0.6 was used from the spectral index map between A6 cm 
and A22cm (Fig. 13 I, which is less affected in this region by 



ionisation losses. The corresponding fluxes and intensities were 
used from the A22 cm map. 

The pixel based technique of calculating the magnetic field 
strength is only applicable for the total intensity emission as an 
estimate of the total magnetic field strength. The polarisation and 
its implication to the ordered magnetic field in M82 will be dis- 
cussed in Paper II. The number density ratio of protons to elec- 
trons was assumed to be constant at a value of Kq = 100. This 
assumption may not be fulfilled due to the possible high frac- 
tion of protons suffering less energy losses due to their higher 
mass than the electrons on their path away from the disk. This 
would lead to Kq > 100 and therefore the calculated magnetic 
field strength gives only a lower limit. The uncertainty here only 
scales approximately with the fourth root of the assumed ratio, 
so that even for values between Kq = 100 to Kq = 1000 the error 
is less than 30%. 

The map resulting from these calculations is presented in 

Fig. 
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Averaging the magnetic field strength over the whole 
galaxy results in a mean equipartition field strength of B eq * 
35 \iG for the assumed geometry. This is high but still close to 
the values for spiral galaxies, which were calculated as having 
9 ± 3 \iG from a sample of 74 galaxies (Ni klas et al.| 1995) For 
galaxies with higher star formation rates like M5 1 the average 
field strengths are 20 \iG ( Fletc her~et al.|[20TT| Looking into 
the central region of M82 gives an even higher magnetic field 
strength of B eq « 98 |iG This may be uncertain by an order of 
magnitude for the above mentioned reasons but is of the same 
order as the value derived by Klein et al. ( 1988) (B eq « 50 ioG). 
Lacki & Beck (in prep.) suggested a revised equipartition for- 
mula for starburst galaxies, which includes pion decay processes 
and the generation of secondary electrons. They calculated a 
magnetic field strength of B eq = 250 iiG, but this value was de- 
rived by using the integrated flux over the whole galaxy and not 
distinguishing between the core and halo emission as well as us- 
ing a constant scalelength of / = 500 pc and is therefore higher. 

The calculated average halo field strength for M82 is B eq « 
24 \iG. Some rising gradient is visible towards the edges of the 
map. Since the resulting field strength is proportional to (j) , the 
assumed geometry influences the values. A more spherical ge- 
ometry would result in a more homogeneous magnetic field gra- 
dient over the halo. Since the geometry of the synchrotron emit- 
ting medium is mostly unknown, the errors are quite huge. But 
anyhow it has to be stressed that the cylindrical geometry seems 
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Fig. 16. Maps of the average equipartition strength of the mag- 
netic field for the assumed cylindrical geometry of the syn- 
chrotron emitting medium using the intensities at A22cm and 
the spectral index in the halo from Fig. T31 The overlaid are the 



soft X-ray (0.2-1 keV) contours from available XMM-Newton 
archive data (Wezgowiec et al. in prep.). The resolution of the 
image is 15". 



to be the preferred one for M82 due to its correlation to other 
wavelengths like Ha (Ohyama et al. 2002) or X-ray (Strickland 
|& He ckman 2009) These Ha observations also showed two dif- 
ferent components of the outflow cone, one with a small open- 
ing angle and a cylindrical geometry and the other with a larger 
opening angle and a plume-like geometry. The latter one would 
reduce the magnetic field strength in the mentioned regions. 



6. Cosmic ray electron losses 

Cosmic ray electrons, which are assumed to be produced in 
star-forming regions mainly suffer energy losses from five dif- 
ferent processes: Synchrotron radiation, the inverse Compton 
(IC) effect, non-thermal bremsstrahlung, ionisation, and adia- 
batic losses. In order to get an estimate of their contribution 
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to the total losses, the different timescales were calculated for 
the derived values for the core and the halo region. This can be 
achieved using the equations given by |Longair (2010) All de- 
rived values are listed in Tab. [7] 

Electrons gyrating in a magnetic field B(\xG) of a given 
strength and radiating at a frequency v(Hz) have the energy 
^(GeV) according to the equation: 



E « 2.5 



urMir 



The synchrotron lifetime t sy „ can then be calculated using: 



,35 • lO^-'B^yrs, 



(15) 



(16) 



To calculate the inverse Compton lifetime t ic an energy density 
Uph (erg cm 4 j is needed: 



Tjc * 3.55 ■l(T A E- x U ph - x yrs, 



(17) 



This can be calculated assuming black body radiation with the 
Stefan-Boltzmann law 



U 



ph 



(18) 



with cr being the Stefan-Boltzmann constant and T the dust 
temperature.. Using a value of T c = 48 K for the centre from 
Kaneda et al. ( 20 1 Q)> leads to an energ y density of U p u, 

A H1 1 . 1 f»-8 am D^vnccal at al I Of>1 f\ I frainH iralnoc Viohi 



4.01310 ergcm .|Roussel et al. ( 2010 1 found values between 



12 K and 50 K for the halo regions. Using an average for the halo 
of 7), = 30 K, which was also estimated by Krugel et al. ( 1990) 



brems 
3\ 



leads to an energy density of U p h,h = 6.124 • 10~ y ergcm"-'. 

To estimate the non-thermal bremsstrahlung lifetime t;, 
and the ionisation lifetime t, „, the ISM particle density n (cm 
is needed: 

Tbrems * 4 ■ lOV 1 yrs (19) 



yrs 

9.5 ■ 10V'£yrs 



Using n c = 250 cm 3 for the core region from 



WeiB et al. 



(20) 
( |2001| ) 



gives a lower limit since |Naylor et aL ( 2010| ) reported particle 
densities of up to several 10 4 cm* 3 . This would reduce Tbrems an d 
Ti on drastically making them the dominant loss processes for the 
core region. For the halo the measured HI mass of 0.75 ■ 10 9 M Q 



from Yun et al. ( 1994| ) is used. Assuming a spherical volume 
with a radius of 5 kpc leads to a particle density in the halo of 
n h = 0.058 cm- 3 . 

The adiabatic losses are only dependent on the derivative of 
the velocity in direction of the minor axis: 



Jdv 



(21) 



The values i n the literature diffe r between velocities of v = 
^OOkms" 1 ( |Reuter et al.||l992t ! for a non reaccelerated cos- 



( Seaq uist & Odegard 



mic ray transport and v = 4000 km s 

1991) for one with reacceleration, but with the close correla 
tion between the radio continuum morphology and the kinematic 
studies using charged particles, we can assume that the mag- 
netic field is frozen into the ionised medium and moving at the 
same speed. Therefore we can take values of v = 600 km s -1 at a 
height above the disk of z = 450 pc from |Bland & Tully| fl988 ) 
w hich is very clo se to the value derived for starbursting galaxies 
by Martin ( 1999 ) of v = 500 km s _1 . The cosmic ray bulk speed 
is an addition of the outflow speed of the charged particles and 



Table 7. Overview of the timescales of the cosmic ray electrons. 
All values were calculated for the A22cm observations. The mid- 
dle column represents the loss timescales for the core region and 
the right column the ones for the halo. 





Bcore 


98 uG 


Bhalo - 


24 uG 


£(GeV) 


0.94 


1.90 


T™ (yrs) 


9.21 


10 5 


7.62 


10 6 


Tic (yrs) 


9.40 


10 3 


3.05 


10 4 


Tbrems (V~S) 


1.60 


10 5 


6.90 


10 s 


Tion (yrs) 


3.57 


10 4 


3.12 


10 8 


Tad (yrs) 


2.20 


10 6 


2.20 


10 6 


Tcool (yrs) 


9.30 


10 3 


3.04 


10 4 


Tesc (Yrs) 


1.17 


10 5 


1.12 


10 s 


Tn (yrs) 


2.00 


10 5 


8.62 


10 8 



the Alfven speed v^. The latter one can be calculated using the 
standard textbook formula 



v A = 2.18w f , 5 fikms 1 . (22) 
This results in a value of Va = 10 km s using n e = 30 cm 



-3 



from McKeith et al. ( 1995 i showing that in this environment the 
Alfven speed is contributing negligibly to the bulk speed. 

The adiabatic losses are still 1-2 orders lower than the inverse 
Compton and ionisation losses which are dominating the loss 
processes in the core as well as in the halo. In addition to the al- 
ready high inverse Compton losses in the core region, ionisation 
losses are roughly of the same order supporting our interpreta- 
tion of the lower fluxes in this region at longer wavelengths. We 
can now calculate if the cosmic rays can escape the galaxy using 
the cooling timescale r c00 /, which can be estimated for a contin- 
uous injection of cosmic rays by combining the synchrotron and 
the inverse Compton timescale: 



-l 

' cool 



-1 

' syn 



-\ 
l lC 



(23) 



This timescale can then be compared to the escape time for the 
cosmic rays due to advection with a galactic wind assuming a 
constant acceleration of the galactic wind from the core region 
into the halo: 



1.02 /zv-'yrs, 



(24) 



3.3.2 1, and v 



where h is the scaleheight of the disk, which was calculated to 
h = 70 pc from the average of the first component of the two 
component fit to the scaleheights for the core region and to h - 
670 pc from the exponential fit to the halo (Sect, 
the bulk speed of the outflowing cosmic ray gas. 

Comparing now the different timescales shows that syn- 
chrotron losses are not the dominant loss process for the cosmic 
ray electrons in the core region even assuming energy equipar- 
tition. Since t coo \ is an order of magnitude smaller than T esc , the 
cosmic ray electrons cannot escape the recent starburst directly 
from the star bursting region into the halo within a normal galac- 
tic wind. Even more important is that the inverse Compton losses 
in the halo, which are dominating the escape timescale, are still 
so high that the cosmic rays are not able to escape the galaxy. 

This shows that the existence of the large radio contin- 
uum halo around M82 is lacking an easy explanation and more 
complicated processes like secondary electrons originating from 
pion decay and multiple star-bursting episodes have to be taken 
into account. 



The IVERITAS Collaboration et al.| ((2009) detected y-ray 
emission from M82 with a luminosity of L\c e v ~ 1.9 ■ 
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10 40 ergs _I at 1 GeV. Scaling this value down to lOOMeV us- 
ing their spectral index of a ~ 2.2 results in LioomcV ~ 8.3 ■ 
10 39 ergs _1 . This value is of the same order as the predicted 
value of LiQQMev = 3 ■ 10 39 erg s~' using L mMeV = 2 ■ 10~ 5 L IR 
( Thompson et al.|20 06) in case of pion decay. This would pro- 



duce secondary electrons in the star-forming regions with the 
same spectrum as the protons and additionally flattening the 
spectrum towards lower frequencies. From the particle den- 
sity n one can calculate the pion loss timescale (Mannheim & 
|Schlickeiser| 1994[> : 



2- 10 : 



(250 cm- 3 ) 



yrs 



(25) 



If this is shorter than r esc the starburst region is a proton 
calorimeter and pion decay has to be taken into account. This 
is the case for the calculated values. But care has to be taken 
here with the estimated values for n since no filling factor was 



included in Eq. 25 It is very unlikely that these secondary elec- 



trons are uniformly produced in the whole core region due to the 
high dumpiness of the neutral as well as ionised gas. 



6.1. Building up the halo 
Using a simple hydrostatic 



Thompson et al. 



region can be ca 



equilibrium approach from 
(2006) the magnetic field strength for the core 
culated with 



J HE 



8tt 2 G. 



v^ 1 



(26) 



where G is the gravitational constant, rj the ratio in which the 
gravitational and magnetic energy are in equilibrium, f g the 
complete gas fraction, and "E g the surface density. Assuming a 
value of 77 = 1 to set an upper limit for the magnetic field 



strength, a gas fraction of f g = I (Forster Schreiber et al. 



and a surface density of 0.85 gem (Kennicutt||1998 1 results 



2001 



in an average magnetic field strength of Bhe = 1.13 mG. This 
would reduce t jv „ to 2.37 • 10 4 yrs and Tjc to 3.19 • 10 4 yrs mak- 
ing synchrotron losses the dominant loss process in this region. 
Using these values and an average bulk speed for the cosmic rays 
of v c = 1000 km s in the core region results in a travel distance 
of s c = 14 pc within the loss timescale, which is too short to 
build up the observed extended synchrotron halo. 

A possible explanation results from a scenario with a filling 
factor of the magnetic and radiation field similar to the ionised 
HII medium of f e = 1.9%. Using our value of B eq = 98 \iG for 
the core region results in an average magnetic field strength of 
B x 5 mG for the individual HH-regions and supernova remnants 
yielding a synchrotron loss timescale of T syn = 2.43 ■ 10 3 yrs. 
Using the same value for the bulk speed as above gives a travel 
distance of s c = 2.5 pc. [Muxlow et al. (1994) found an average 
size of 2.4+ 1 .0 pc for the supernova remnants in M82. This value 
is by more than one order of magnitude smaller than the one for 
supernova remnants in the Large Magellanic Cloud or the Milky 
Way and therefore magnetic field strength in the mG regime are 
a reasonable value. It was interpreted as being smaller due to its 
young population, but it can also be smaller due to the higher 
pressure of the medium surrounding those remnants and con- 
straining their expansion. This would additionally support the 
arguments for a low filling factor of the ionised medium and the 
magnetic and radiation field. 

Fig. 17 shows the dependency of the cosmic ray energy vs. 
their travel times in case of an exploding supernova inside an 




10 9 

Frequency (GHz) 



Fig. 17. Plot of the travel distance for the cosmic ray electrons 
for the combination of synchrotron and inverse Compton losses 
vs. their emitted frequency inside the core region under the as- 
sumption of a magnetic field of B « 5 mG and a bulk speed 
of v c = 1000 km s . The straight line shows the average size 
for the compact sources from |Muxlow et a l. (1994) with their 
errors marked as shaded. The two lines cross each other at 
v = 1.38 GHz. 



HH-region using the above mentioned field strength and veloc- 
ity at frequencies between 100 MHz and 10 GHz. A comparison 
with the average compact source size (dashed line) shows that 
only electrons radiating at frequencies v < 1.38 GHz can es- 
cape the HH-regions. Since synchrotron emission from the halo 
is only visible at A22 cm and A92 cm this explanation fits the ob- 
servations well. 

The pressure to keep the HH-regions and supernova rem- 
nants confined would then originate from the hot X-ray gas as 
well as the molecular and neutral gas filling the cavity in be- 
tween them with a filling factor close to 100%. The gas density 
and photon pressure is up to 10 4 times higher than those in the 
Milky Way on larger scales resulting in a much higher filling 
factor for these gas phases. A simple estimate using pressure 
equilibrium between the ionised gas in the supernova- and HII- 
regions with a temperature of 10 5 K and the X-ray gas with an 
assumed temperature of 10 7 - 10 s K leads to an energy density of 
50 -500 eVcm~ 3 . This is on the same order as the value derived 
by |Klein et al. 



(1988) (eOeVcirT 3 ), but at least by a factor of 



100 lower than the energy density of the photon field. Including 
the high molecular and neutral gas densities of up to 10 4 cm~ 3 
results in pressures of the same order as the photon field. To ver- 
ify such a distribution further one would need observations with 
an improved resolution at X-ray and mm wavelengths, which are 
not available with current instruments. 

Since we showed that electrons produced by supernovae in 
contrast to their accompanying protons already lose a large part 
of their injected energy over the first parsecs of their travel, the 
ratio between the energy densities in protons and electrons in- 
creases towards the halo. Additionally, the hot Ha gas filling the 
whole galaxy is charged and will couple the gyrating electrons 
to its kinematics together with the magnetic field as soon as they 
leave the high magnetic field in the HH-regions and inside the 
supernova shells. This explains the strong correlation between 
the Ha-, PAH + -, and radio continuum morphologies and still 
account for the observed indications concerning pion decay and 
secondary electrons. 
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Using a simple estimate for the travel distance in the halo 
with a cosmic ray bulk speed of v = 600kms~ 1 and the inverse 
Compton losses as the dominant loss process results in a travel 
distance of s = 19 pc. This is still by at least two orders of mag- 
nitude lower than the observed kpc extent of the synchrotron 
halo. A decrease of the dust temperature in the halo directing 
outwards was observed by Rous sel et al.| ( [20T0| , which would 
lower the calculated values for the radiation field and would in- 
crease the travel distance by decreasing the inverse Compton 
losses. Since the radiation field is proportional to T A and there- 
fore directly proportional to the travel distance, the dependence 
of the energy losses on the assumed dust temperature is signif- 
icant. Using T/, = 10 K from |Roussel et alT| ( |20 1 0) > as a lower 
limit for the dust temperature, results in a photon field strength 
of U p h,h - 7.56 • 10~ n ergcirT 3 and therefore a travel distances 
of s ~ 1.5 kpc. This shows, that arguing for a continuous injec- 
tion of cosmic rays from the centre into the halo would need dust 
temperatures of T < 10 K, which does not fit the observations. 
A scenario, which does not interfere with those, is the build-up 
of the halo within multiple starbursting periods, where the radia- 
tion field varies strongly and allows the fuelling of the halo with 
fresh material. 



6.2. Multiple starbursting periods 

Forster Schrei ber et al.| (|2003 ) describe a scenario of at least two 
successive starburst epochs for M82 about 10 and 5 Myrs ago in 
addition to the recent one, where the material was released into 
the halo and is accelerated by shocks and the high pressure from 
the core region ( |Kaneda et al.||2010| This would also explain 
the observed kpc extent of the synchrotron halo. Within periods 
where the radiation field due to less star-formation declines, the 
cosmic rays can leave the inner region of the galaxy and build 
up a halo of the observed size. Since tjq is proportional to U p h, 
this would need variations of the radiation field over an order of 
magnitude over the time of a starburst period to reach heights 
above the disk of several kpcs. |Mas-He sse (1992 ) showed that 
variations in the far-infrared luminosity of this order are reason- 
able for young, compact starbursts like M82. 

Such an evolutionary scenario would not only explain the 
extent of the radio halo, but also the abundance of PAH + in the 
halo, which would if expelled into the halo from the current star- 
burst have been destroyed due to the strong radiation from the 
photon field ( |Bendo et~a l. 2008) Furthermore, the hot gas cap 
visible far to the north in Fig. 16 might be a remnant of an older 



On the other hand the recent high-resolution Hot-observations 
show small shocks inside the outflow cone suggesting that at 
least on these scales the cosmic ray electrons reside inside a 
small-scale turbulent field. This cannot rule out that on larger 
scales the outflow may still be polarised. |Lerche & Schlickeiser 



starburst phase at the interface of the intergalactic medium. 

The difference between the extents of the synchrotron halo at 
A22/A92 cm and A3/A6 cm is now explainable. The recently in- 
jected higher energetic cosmic rays would not reach the halo due 
to their high synchrotron and inverse Compton losses in the com- 
pact HII- and supernova-regions in the core. This would even 
be the case for an episode between two starburst periods, since 
the pressure of the photons only delivers ~ 10% to the pressure 
balance in the core between the compact regions and their sur- 
roundings. In this case the cosmic rays could be expelled into the 
halo due to their increased inverse Compton timescale, which is 
not possible during a starburst period like M82 is encountering 
now. 

Possibilities for particle reacceleration and diffusion were 
discussed by Seaquist & Odegard (1991) While reacceleration 
by small-scale turbulence would lead to depolarisation inside the 
beam, large scale shocks originating from an epoch of starburst 
activity can still reaccelerate the electrons in the halo and lead to 
observable polarisation degrees as seen by Reuter et al. (1994) 



( 1980) suggested a diffusion dominated model which would re- 
sult in a flat spectral index for the halo. This would not interfere 
with the observations and therefore may be the preferred one. 
Since the analysis of the polarisation data is going to be dis- 
cussed in Paper II, we refer to this for a detailed discussion. 

7. Summary 

A collection of datasets from the VLA at A3 cm and A6 cm and 
from the WSRT at A22 cm and A92 cm was reduced to carry out a 
multifrequency analysis of the starburst galaxy M82. The phys- 
ical parameters in the core region as well as its influence on the 
halo were investigated and connected to the evolution scenario 
of the galaxy. 

Advanced data calibration techniques for radio continuum 
data are needed to achieve dynamic ranges sufficient for the sci- 
entific analysis of strong extended sources like M82. With the 
newly proposed technique in this publication we were able to 
reach dynamic ranges of ~ 50000 for A22 cm observations and 
even ~ 10000 for A92 cm. 

These high dynamic range images allowed us to look into 
the morphology of the halo at these long wavelengths. While 
the extent of the halo at A3 cm and A6 cm is limited to the inner 
1 kpc radius from the core, diffuse emission up to 4 kpc radius 
shows up at A22 cm and A92 cm. A large difference in the extent 
of the northern and southern outflow is recognised. While the 
northern emission reaches up to 4-5 kpc into the halo, the south- 
ern one can only be traced up to 2.5 kpc. This asymmetry is also 
visible in the X-ray observations (Fig. 16 1. A possible explana- 
tion would be a movement of M82 towards the group centre in 
the south causing ram pressure effects like in clusters of galax- 
ies, but such an evolution scenario contradicts the results by |Yun| 
|etaLldl993a) . 

The northern outflow cone shows a lack of emission at long 
wavelengths in contrast to the two spurs surrounding it and steep 
spectral indices of up to -2.5 in its central part. This indicates 
that the dominant loss processes for the outflow cone in the north 
are either adiabatic or synchrotron losses. Numerical studies by 
Miyake et al. (2007 ) confirm this in the case of NGC253, where a 
similar but not as well pronounced morphology was recognized 
by |HeesenetaLlp009) . 

A fitting to the scaleheights lead to very different results than 
in spiral galaxies. While for spiral galaxies typically a thin and 
a thick disk is found with scaleheights of 0.3 kpc and 1.8 kpc 
respectively (Krau se|201 1 ij the values for M82 were lower by a 
factor of ~ 3. Additionally an overall small gradient for the spec- 
tral index was found for the halo between A22 cm and A92 cm 
with a = -1.2 1.5, while an increase is expected for galax- 
ies where inverse Compton and synchrotron losses dominate the 
cosmic ray losses. Using an energy equipartition ansatz between 
the magnetic field energy and the cosmic ray particles results in 
an average field strength of 24 for the halo and 98 |xG for the 
core region. 

An integration over the emission of the core and the halo 
showed that free-free absorption is lowering the flux of the core 
at wavelengths from A22 cm on, while the halo emission is com- 
pletely free of this feature at least up to A92 cm indicating a high 
contribution of ionisation losses in the core region. From a fit 
to the data of the core an average electron density of 22 cm" 3 
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could be obtained leading to a filling factor of ~ 1.9% for the 
ionised hydrogen. These fits are limited to just five datapoints 
and therefore have to be taken with care, but additional obser- 
vations between A22cm and A92cm with the WSRT and even 
at lower wavelengths with LOFAR will give better constrains on 
this. 

A calculation of the loss timescales shows that at the current 
stage the inverse Compton and ionisation losses are dominating 
the expansion of the cosmic rays in the core region confirming 
the previous statement. Additionally the recent y-ray detections 
support the thesis of M82 encountering pion decay in the core 
region indicating very high particle densities, but lack the reso- 
lution to determine the exact origin of these particles within the 
galaxy. 

To account for the above observations we pointed out a sce- 
nario where HH-regions and supernova remnants are dominat- 
ing the radio flux in the core region. Cosmic rays are generated 
inside these small regions, where the medium is so dense that 
secondary electron production via pion decay is possible and 
cosmic ray electrons are suffering high losses due to magnetic 
fields. Using the same filling factor for the magnetic field as 
for the ionised medium, the magnetic field strength inside these 
compact regions was estimated to 5.16mG. This makes syn- 
chrotron losses the dominant loss process inside these regions 
and the travel times for these particles at A < 6 cm are smaller 
than size of these regions. This explains the difference of the ex- 
tent of the synchrotron halo between frequencies of A3 cm/A6 cm 
and A22 cm/A92 cm. Higher energetic particles lose their energy 
faster than lower energetic ones and are not propagated out into 
the halo. Since the pion decay mostly occurs inside these dense 
regions, it can not help to produce electrons less affected by these 
high synchrotron losses. 

In contrast to this heavier charged particles like protons or 
molecules are less affected by the magnetic field and therefore 
can escape these small regions more easily. The pressure to con- 
strain these regions to sizes of several pc can be generated by 
the hard X-ray emission and the neutral hydrogen and molecu- 
lar gas filling the whole core region of the galaxy. As soon as 
charged particles reach the lower density regions surrounding 
the star-forming regions they become coupled to the kinematic 
flow of the lower density X-ray gas and may leave the core re- 
gion, which explains the similarity in the morphologies of Ha, 
PAH + and radio. 

The inverse Compton losses are absolutely dominating the 
loss processes in the halo. Using a simple scenario, where cos- 
mic rays are propagated from the core region into the halo results 
in travel times of only 19 pc. To argue for a constant injection of 
cosmic rays into the halo, dust temperatures below 10 K would 
be needed, which do not fit the observations. This gives evidence 
for a scenario where M82 ran through several starburst periods in 
the last 5-10 Myrs, where the photon field decreased by one or- 
der of magnitude and it was possible to eject cosmic rays into the 
halo. The small spectral index gradient in the halo indicates that 
different propagation effects than in normal spiral galaxies are 
dominant. It cannot explain a halo where a continuous propaga- 
tion of cosmic rays suffering from inverse Compton losses is fu- 
elling the halo. Instead we propose that the synchrotron halo was 
produced by multiple ejections of cosmic rays and the observed 
spectrum is a superposition of several populations of those. 

8. Conclusions 

The described scenario could imply that the magnetisation of the 
intergalactic medium by starbursts is limited to intervals shorter 



than the actual time of the starburst since losses for the cosmic 
rays are too high to let cosmic ray electrons escape the starburst- 
ing galaxy as soon as a strong photon field and high densities 
have been reached. Only heavier charged particles might still es- 
cape the galaxy due to their longer loss timescales. This would 
have a significant influence on the proton to electron ratio in 
the intergalactic medium and significantly influence theoretical 
predictions of the propagation of cosmic rays and the strength 
of magnetic fields in the intergalactic space. For the first time 



Kotarba et al. (201 1 1 showed in simulations of three interacting 
galaxies, that the whole group medium can become magnetised 
over a timescale of several hundred Myrs. 

This evolution scenario and the discovered high dumpiness 
of the medium in M82 even has influence on predictions for the 
early Universe. Due to the overall increased star-formation rate 
and number of starburst galaxies the release of magnetic fields 
into the intergalactic medium may not be as easy as predicted 
by current simulations. Additionally the high pressure environ- 
ments in such galaxies might sufficiently suppress the expansion 
of supernova shells. Hence the dominant mechanism for mate- 
rial leaving the originating region are large scale shocks from 
only several supernova explosions inside an HH-region. To ver- 
ify our evolution scenario further numerical MHD-simulations 
are needed using the discussed loss processes together with the 
constrains given here and in several other publications. 

With LOFAR becoming more and more scientifically usable 
and complementary observations from the WSRT and VLA with 
wavelength of up to A92cm for a lot of nearby galaxies the 
free-free absorption at long wavelengths is a viable technique 
to estimate the thermal electron density and the filling factor 
of the ionised medium of star-forming spiral and starbursting 
galaxies. The integrated spectrum over whole galaxies does not 
sufficiently account for the difference between the complicated 
physics in the star-forming regions and the halo and interarm ar- 
eas. The mentioned observations are able due to their good res- 
olution to differentiate between these regions and provide new 
input for theoretical predictions and numerical simulations. 
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